Plant Biotechnology Journal

The transition from hulled to free-threshing grain was a pivotal event in wheat domestication, enabling efficient harvesting and processing. Threshability in tetraploid wheat is controlled primarily by the Q locus and two Tenacious glume (Tg) loci on chromosomes 2A and 2B, yet the molecular basis of the major Tg1-B locus remains incompletely characterized. Here, we phenotyped a durum wheat x wild emmer wheat (WEW) recombinant inbred line (RIL) population across two field environments and performed QTL analysis for glume tenacity (TG), threshability ratio (THRR), and seed number per spike (SDNPS). A total of 19 significant QTLs were detected across six chromosomes. The largest-effect loci for both TG and THRR co-localized on chromosome 2B, with LOD scores up to 14.22 and phenotypic variance explained up to 31.2%, corresponding to the previously described Tg1-B locus. To validate this QTL, the donor RIL was backcrossed three times to Svevo to generate a near-isogenic line, NIL-65 (BC3F5), confirmed by whole-genome skim sequencing to carry a homozygous WEW introgression at Tg1-B. A segregating BC4F2 population derived from NIL-65 confirmed that plants homozygous for the dominant Tg1-B allele displayed significantly higher glume tenacity and intact glume morphology compared to tg1-B sister lines, which exhibited basal glume cracking characteristic of the free-threshing phenotype. Genotyping-by-sequencing delimited the causal interval to an approximately 11 Mb introgression on chromosome 2B. These results confirm the major role of Tg1-B in determining glume tenacity in tetraploid wheat, provide a validated near-isogenic germplasm resource, and lay the foundation for fine-mapping and functional characterization of the underlying gene(s).

Fitness costs of plant disease defence are often subtle and difficult to quantify. In this study, we therefore used comparative high-throughput phenotyping in two independent facilities to assess growth, morphology and physiology of potato (cv. Desiree) with high time-resolution monitoring different defence mechanisms under pathogen-free conditions. Plants were either treated weekly with the resistance inducers {beta}-aminobutyric acid (BABA; 10 mM) or potassium phosphite (KPhi; 36 mM) or comprised six transgenic lines expressing late blight resistance genes (single Rpi genes or a three-gene stack) or reduced jasmonate perception (StCOI1-RNAi). Over four weeks, image-derived traits revealed consistent cross-facility effects for plant height and colour: BABA treatment increased plant height but reduced canopy area and induced a paler greenness signature, whereas KPhi caused minimal and transient growth effects. Chlorophyll fluorescence at the NaPPI facility indicated reduced vitality (Rfd_Lss) in BABA-treated plants and increased Rfd_Lss following KPhi, while maximum PSII efficiency was largely unchanged. Several transgenic lines showed somewhat reduced above-ground biomass. Enzyme activity profiling produced distinct treatment and genotype signatures, but was strongly modulated by facility conditions that overrode these specificities. Overall, high-throughput phenotyping robustly detected subtle growth-defence trade-offs across platforms. HighlightHigh-throughput optical phenotyping validated across two independent research facilities reveals that stacked resistance genes and resistance inducers in potato trigger subtle growth trade-offs. Graphical abstracts O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=97 SRC="FIGDIR/small/713143v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@89df47org.highwire.dtl.DTLVardef@1a1ce64org.highwire.dtl.DTLVardef@1f52f0dorg.highwire.dtl.DTLVardef@1e41c35_HPS_FORMAT_FIGEXP M_FIG C_FIG Experimental timeline for high-throughput plant phenotyping platforms. Created in BioRender. Poque, S. (2026) https://BioRender.com/nmkve7g

The narrow genetic base of cultivated tomato (Solanum lycopersicum L.) represents a major constraint on crop improvement, necessitating the exploitation of wild relatives to broaden allelic diversity. Here we present SABER (Solanum lycopersicum Allele Biodiversity Enriched Resources), a novel eight-founder Multiparent Advanced Generation Intercross (MAGIC) population that, for the first time, incorporates the Galapagos wild relative Solanum cheesmaniae as a founder alongside seven elite S. lycopersicum lines. Following a structured crossing scheme and Single Seed Descent advancement, F6 recombinant inbred lines were genotyped at 5,850 high-confidence SNP markers using Single Primer Enrichment Technology (SPET). Population structure analyses confirmed low residual heterozygosity, limited substructure among offspring, and successful introgression of S. cheesmaniae alleles across all twelve chromosomes. Mapping performance was validated through three Mendelian traits with known genetic determinants, all of which resolved to genomic positions consistent with the literature. QTL mapping for quantitative agronomic traits identified known loci for fruit epicarp and flesh color, and two novel QTL for days to flowering, number of leaves before flowering, and soluble solids content. Together, these results demonstrate that SABER is a powerful and reliable platform for high-resolution QTL mapping and candidate gene discovery, and establish a replicable framework for integrating wild germplasm into multiparental tomato breeding resources

Course-based Undergraduate Research Experiences (CUREs) have emerged as a transformative approach to science education, expanding access to authentic research opportunities beyond the traditional undergraduate research assistant (URA) training. By embedding research into a curriculum, CUREs engage a broad and diverse population of students in a classroom environment that emphasizes experimental design, data analysis, and scientific communication. However, this has been difficult to develop for fields such as plant synthetic biology due to the long timescales of plant transformation. One avenue around this problem is to utilize a recent innovation that enables high throughput and rapid screening of gRNA efficacy by leveraging viral-based delivery of guide RNAs (gRNAs). In this work, we develop and validate a CURE with undergraduate students at Colorado State University (CSU). Students worked in teams to design and test efficacy of gRNAs targeting a Cas9-based transcriptional repressor to different regions of the promoters of the three GIBBERELLIN INSENSITIVE 1 genes (GID1a, GID1b, and GID1c) in Arabidopsis thaliana. Over the semester, students generated and analyzed gene expression data to understand the efficiency of twelve new gRNAs. We further validated CURE student-identified gRNAs with an undergraduate research assistant (URA) that assessed target gene expression and phenotypic outcomes in stable transgenic lines expressing SynTF constructs with the strongest gRNAs from the class. We further describe the curriculum structure to facilitate adoption at other institutions and present student-generated datasets demonstrating the utility of ViN-based screening for identifying effective SynTF gRNAs for plant functional genomics and engineering. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=111 SRC="FIGDIR/small/715601v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@13869f5org.highwire.dtl.DTLVardef@b469feorg.highwire.dtl.DTLVardef@9aa51borg.highwire.dtl.DTLVardef@cdc129_HPS_FORMAT_FIGEXP M_FIG C_FIG

Sacha inchi (Plukenetia volubilis L.) is an emerging woody oilseed crop prized for its high alpha-linolenic acid (ALA) content. Despite its nutritional and economic value, the lack of high-quality genomic resources has hindered genetic improvement and the elucidation of its unique polyunsaturated fatty acid and lipid biosynthetic pathways. In this study, we report a high-quality, chromosome-scale genome assembly of sacha inchi with a total length of 710.62 Mb, integrated from Illumina, PacBio, and chromosome conformation capture (Hi-C) technology. The genome harbors 37,570 protein-coding genes, and 379.86 Mb (53.45%) of repetitive sequences. Phylogenomic analysis reveals that sacha inchi diverged from its closest relative Ricinus communis, [~] approximately 36.2 million years ago. Comparative genomics indicates that sacha inchi experienced only ancient whole genome duplication events. To elucidate the mechanisms governing ALA biosynthesis and triacylglycerol (TAG) accumulation in sacha inchi seeds, we performed temporal transcriptome profiling across six seed development stages. Our findings demonstrate that high TAG content is primarily driven by the sustained expression of biosynthetic genes and low activity of degradation genes during mid-to-late seed development. Notably, while genes encoding stearoyl-ACP desaturases (SADs) maintain the precursor pool, the expression of genes encoding fatty-acid desaturase 2 (FAD2) and fatty-acid desaturase 3 (FAD3) is positively correlated with the final accumulation of C18:2 and C18:3 fatty acids. We also identified lncRNAs as potential epigenetic regulators of these key pathways. This high-quality genome provides a critical foundation for elucidating the molecular mechanisms of seed growth and development in sacha inchi.

Concerns over climate change have intensified the demand for stress resistant crops like hybrid wheatgrass (HWG; Elymus hoffmannii, StStStStHH), a perennial forage species known for its exceptional salt and drought tolerance. However, hexaploidy and high heterozygosity have complicated efforts to resolve its genomic structure and evolutionary history. Here, we present high-quality, haplotype-resolved, chromosome-level genome assemblies for HWG (CDC Saltking) and its putative progenitor, bluebunch wheatgrass (Pseudoroegneria spicata). By integrating PacBio HiFi and ultra-long Oxford Nanopore sequencing with Hi-C scaffolding, we assembled the 10.7 Gb HWG genome into 21 pseudochromosomes per haplotype. Our phylogenomic analysis redefines the origin of the H subgenome, positioning it as an intermediate between Old-World Hordeum marinum (sea barley) and Hordeum brevisubulatum. Notably, we identified significant chromosomal rearrangements, including a unique duplication on St chromosome 4. Transcriptome analysis across multiple tissues revealed a pronounced expression dominance of the H subgenome. This dominance was not associated with reduced LTR density, suggesting that selective pressures for rapid adaptation of the latest subgenome entrant may drive its dominance. Finally, using the f-branch statistic, population genomic analysis of 189 accessions representing eight Elymus and Pseudoroegneria species revealed extensive reticulate evolutionary relationships and identified P. spicata as a major, asymmetric genetic donor within the wheatgrass complex. These resources provide a foundational framework for future genomic research and genetic improvement in grasses and for the introgression of stress-tolerance traits into cereal crops such as wheat. Key MessagesDevelopment of world-first high-quality chromosomal-level haplotype-resolved genome assemblies of hexaploid HWG and diploid progenitor, Pseudoroegneria spicata, enabled the identification of the subgenome origins. This study resolved the evolutionary placement of the St genome and clarified the history of polyploidization and hybridization in HWG. Homeolog expression bias in the H subgenome likely reflects selective pressure favoring greater gene retention and upregulation of functionally important genes, thereby enhancing hybrid fitness. Population structure analysis distinctly differentiates P. spicata, E. repens, E. hoffmannii from other European Pseudoroegneria species. The findings reveal the complex patterns of interspecific gene flow and population dynamics within the Elymus and Pseudoroegneria species.

Common buckwheat (Fagopyrum esculentum Moench) is a globally cultivated pseudocereal with a high nutritional quality and economic value. Due to its self-incompatibility, common buckwheat exhibits a high level of heterozygosity, making genome assembly challenging. Consequently, reference-level haplotype-resolved assemblies of common buckwheat are scarce, hindering research and genomics-assisted breeding. Here, we present a near-complete, chromosome-level, haplotype-resolved assembly of a common buckwheat F1 genotype (named Tuka), generated using a trio-binning approach that integrated parental Illumina short-read data with PacBio HiFi and Hi-C data from Tuka. The Tuka assembly comprises two haplomes, Tuka_h1 and Tuka_h2, both showing high contiguity (contig N50 of 76.68 Mb and 84.57 Mb, respectively), high completeness (assembly sizes of 1.28 Gb and 1.23 Gb with BUSCO scores of 96.9% and 96.8%, respectively), high base-level accuracy (QV of 59.08 and 63.03, respectively), and few gaps (35 and 30, respectively). This near-complete assembly of Tuka serves as a valuable genomic resource for common buckwheat, enabling advanced genomic analyses and accelerating research and breeding using state-of-the-art genomic tools.

Computational tools, including biological language models (LMs), show substantial promise in predicting the impact of genetic variants on plant fitness. However, validating variant effect predictions (VEP) requires experimental populations where genetic variation consists of discrete point mutations rather than segregating recombination blocks. In this study, we generated a novel population of Brachypodium distachyon mutant lines to evaluate the accuracy of VEP at single-base resolution. These lines were advanced through single-seed descent for five generations (M1 to M5), with whole-genome sequencing performed at M2 and M5 and phenotypic measurements recorded at M3 and M4. Using state-of-the-art VEP models, we predicted the functional impact of missense protein-coding variants and gene-proximal non-coding variants. We validated these predictions by estimating the effect of mutations on whole-plant measurements (burden tests) and their probability of fixation from M2 to M5 (purging tests). Among missense variants, the protein LM ESM showed superior predictive accuracy compared to the bioinformatic standard SIFT and the genomic LM PlantCAD. Notably, the relationship between VEP scores and allele fixation suggested a log-linear relationship between VEP scores and variant fitness. Among gene-proximal variants, PlantCAD appeared more accurate than supervised models of regulatory activity, such as chromatin accessibility (a2z) and RNA abundance (PhytoExpr). Collectively, our findings highlight the utility of state-of-the-art VEP tools as predictors of fitness and demonstrate the potential of mutant populations to evaluate computational tools for precision breeding applications.

We present a global soybean haplotype map generated from whole-genome sequencing of 1,278 Glycine max and Glycine soja accessions, comprising 11.37 million SNPs and 2.05 million short insertions and deletions. This map (GmHapMap-II) captures unprecedented worldwide genetic diversity, reflecting the broad extent of the global soybean gene pool. Population structure analyses revealed six geographically distinct subpopulations that affected the linkage and shaped the recombination. The haplotype variation map was used to identify novel genomic regions associated with crude protein content on chromosome 15 that were not detected by a lower SNP density array. LD-based haplotype analysis revealed a superior haplotype for crude protein content. The constructed haplotype map enabled detailed characterization of haplotype diversity and copy number polymorphism at the SCN-associated rhg-1 and Rhg-4 loci, revealing both novel haplotype structures and germplasm lines with elevated CNV relative to previously characterized genotypes. We employed the HapMap matrix for a multi-class variations ML-based genomic prediction approach to predict phenotypes for SCN and catalogued the gene-centric haplotypes in a user-friendly database. The analysis revealed the extent of deleterious alleles present in the soybean germplasm and how breeders have deployed beneficial alleles and purged deleterious ones. The haplotype map will serve as a major genomic resource for trait-based mapping, enhancing efforts in the genomics-enabled development of improved cultivars.

Carmenere is a widely cultivated and internationally recognized grapevine cultivar in Chile, yet genetic variation among its clones remains poorly characterized. Early studies based on SSR and AFLP markers detected limited polymorphism, but these approaches interrogate only a small fraction of the genome, leaving the extent of clonal diversity unresolved. Here, we generated an improved chromosome-scale diploid genome assembly of Carmenere FPS clone 02 and characterized clonal genomic diversity by sequencing 36 biological replicates representing 12 clones maintained in Chile, including heritage selections rescued from old producer vineyards by Vina Santa Carolina as part of its Bloque Herencia conservation program, and commercial nursery-derived clones. Focusing on low-frequency variants and using replicate-aware consensus calling, we identified more than 9,000 private single nucleotide variants (SNVs) and small indels per clone, providing high-resolution markers for clonal identification. Although most variants were located in repetitive or intergenic regions, a subset affected coding sequences, with genes involved in plant-pathogen interactions, transport, and secondary metabolism most frequently impacted. While variant-affected genes associated with wine anthocyanin content, TA, pH, and alcohol percentage were identified, broader phenotypic characterization will be required to assess their biological significance. Overall, this study provides a genome-wide characterization of extant clonal diversity in Carmenere, with implications for clonal selection and genetic resource conservation.

Strigolactones (SLs) are plant hormones regulating shoot branching and symbiotic interactions, but their trace-level abundance limits research and applications. Here, we optimized a Nicotiana benthamiana transient expression system for SL production by tuning agroinfiltration parameters and co-expressing rate-limiting carotenoid biosynthetic genes. Overexpression of Zea mays PSY1 or an Arabidopsis PSY-GGPS11 fusion increased carlactone production over 2-fold and enhanced downstream SL accumulation. Using this platform, we discovered that sorghum cytochrome P450 SbCYP728B35 catalyzes conversion of 5-deoxystrigol to sorgolactone, revealing a previously unknown function. These results establish metabolic engineering of precursor supply as an effective strategy for boosting SL production and demonstrate N. benthamiana as a robust system for pathway elucidation and biotechnological synthesis of bioactive strigolactones.

Pepper (Capsicum annuum L.) is a recalcitrant species regarding shoot regeneration, a trait that serves as a major bottleneck for the application of genetic engineering tools. In this study, comparative genetic analysis between a rare high-regeneration cultivar and a common low-efficiency cultivar identified a single nucleotide polymorphism (SNP) in PHYTOCHROME A SIGNAL TRANSDUCTION 1 (CaPAT1) that determines shoot regeneration efficiency. The T478C SNP in the high-efficiency cultivar converts a stop codon into an Arg codon, leading to translational read-through into the neighboring gene and forming an intact GRAS domain. This SNP-mediated formation of full-length CaPAT1 is essential for its dimerization. Notably, the overexpression of CaPAT1T478C in multiple low-efficiency cultivars, including both hot and bell peppers, significantly improved both shoot regeneration and transformation efficiency in the transformed T0 generation. These findings demonstrate the pivotal role of CaPAT1 in enhancing shoot regeneration and provide a robust strategy to overcome recalcitrance in pepper.

Fast-neutron mutagenesis creates diverse genome-wide mutations, providing a powerful tool for crop functional genomics. Here, we present an expanded genomic and phenotypic analysis of 3,268 fast-neutron (FN)-induced mutant rice lines (Oryza sativa L. cv. Kitaake). All FN lines were whole-genome sequenced, and mutations were identified by alignment in the Nipponbare and KitaakeX reference genomes. We cataloged over 428,000 mutations affecting 78.49% of Nipponbare genes and 70.38% of KitaakeX genes. In silico expression analysis indicates that 575 non-mutated Nipponbare genes are highly expressed and likely essential for viability. Each mutant carries, on average, 68.5 mutations in the Nipponbare alignments or 63.2 mutations for KitaakeX alignments, distributed randomly across all 12 chromosomes with no evident hotspots. FN lines have approximately 8.5% fewer mutations when using the KitaakeX alignment, underscoring the unique contributions of each reference genome and the importance of utilizing both for comprehensive mutation discovery. The majority of mutations are small deletions and single-base substitutions, with deletions predominating in their effect on genes. We found that 74.4% of all transcription factor Nipponbare genes were mutated at least once. Phenotypic characterization of over 2,700 lines revealed a broad spectrum of variation in core agronomic traits (heading date, tiller number, plant height, panicle weight, seed yield components) and other morphological variants of interest. The integration of genomic and phenotypic data through the KitBase platform enabled the identification of candidate genes for several traits of interest. The KitBase website (https://kitbase.ucdavis.edu) has been updated to provide open access to all mutation data and seed stocks, as well as an intuitive query interface, facilitating forward and reverse genetic analyses in rice. This expanded resource enriches the rice functional genomics toolkit and highlights the value of coupling high-density mutation mapping with phenotypic data for rapid gene discovery and crop improvement.

Hi-C data is commonly used for reference-free de novo scaffolding. However, with the rapid increase in high-quality reference genomes, reference-guided workflows are now more practical for assembling large numbers of target genomes without relying on costly and labor-intensive Hi-C sequencing. Recently, a pangenome graph-based haplotype sampling algorithm was introduced to generate personalized graphs for target genomes. Such graphs have strong potential as references for reference-guided contig scaffolding. Here, we present noHiC, a reference-guided scaffolding pipeline supporting key steps of plant contig scaffolding. A distinctive feature of noHiC is the nohic-refpick script, generating a best-fit synthetic reference (synref) from a pangenome graph that is genetically close to the target contigs. This enables the integration of genetic information from many references (up to 48 in our tests) without using them separately during scaffolding. Synrefs showed advantages over highly contiguous conventional references in reducing false contig breaking during reference-based correction. Additionally, nohic-refpick can be combined with fast scaffolders (ntJoin) to rapidly produce highly contiguous assemblies using synrefs derived from pangenome graphs. The noHiC pipeline, used alone or in combination with ntJoin, can generally produce assemblies that are structurally consistent with public Hi-C-based or manually curated genomes. The pipeline is publicly available at https://github.com/andyngh/noHiC. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=82 SRC="FIGDIR/small/712436v1_ufig1.gif" ALT="Figure 1"> View larger version (9K): org.highwire.dtl.DTLVardef@40bd8forg.highwire.dtl.DTLVardef@5d2bbborg.highwire.dtl.DTLVardef@e214a3org.highwire.dtl.DTLVardef@b90b06_HPS_FORMAT_FIGEXP M_FIG C_FIG

T-DNA insertion mutants are widely used to disrupt genes and infer their functions, yet the insertions can also trigger unintended genomic changes that confound phenotypic interpretation. Here, we used T-DNA insertion mutants affecting the major mitochondrial malate dehydrogenase (MDH1) and the heterodimeric NAD-dependent malic enzymes (ME1 and ME2) to examine their functional coordination across photoperiods and irradiance regimes. Under short days, especially at low light intensity, mdh1xme2 mutants were markedly smaller than wild type and, unexpectedly, than the mdh1xme1xme2 triple mutant, and they showed a more pronounced reduction in photosynthetic capacity. ME1 was undetectable in mdh1xme2, implying that the double and triple mutants effectively lack heterodimeric ME and should therefore behave similarly, contrary to what we observed. Whole-genome analysis resolved this discrepancy by revealing that the MDH1 T-DNA insertion in mdh1xme2 is accompanied by a major rearrangement, a 137-kbp duplication downstream of the insertion site, which was absent in the mdh1xme1xme2 triple mutant. This duplication increased gene dosage and elevated transcript abundance across the duplicated interval, while proteomics detected 5 of the 38 encoded proteins, including PEPC1. mdh1xme2 accumulated oxaloacetate-derived amino acids and displayed an altered carbon/nitrogen balance, making PEPC1 a plausible contributor to the exacerbated mdh1xme2 phenotype. Together, our data indicate that a T-DNA-linked structural variant can amplify expression of dozens of genes and intensify phenotypes at specific conditions, thereby affecting the interpretation of genotype-phenotype relationships. Because Agrobacterium-mediated DNA transfer also underpins many genome-editing workflows, our findings argue that structural validation around insertion/editing loci should be considered essential when interpreting T-DNA-derived plant lines.

Fusarium oxysporum FO12 was originally isolated from cork oak (Quercus suber L.) and has been characterised as a highly effective biological control agent of wilt diseases on different crops. FO12 endophytically colonises roots and basal stems of plants, reducing the establishment of the soil-borne pathogen Verticillium dahliae and triggering plant defence-related genes. Here, we report a chromosome-level genome assembly of FO12 using Nanopore and Hi-C data. The 57.60 Mb assembly comprises 14 chromosome-scale scaffolds with centromeres resolved and telomeric repeats detected at 4 of 28 chromosome ends. This high-quality reference genome provides a valuable resource for further research into the use of FO12 in agriculture as a biocontrol agent.

Grain mold severely constrains sorghum [Sorghum bicolor (L.) Moench] productivity and grain quality in subhumid environments. Photoperiod-sensitive flowering plays a key role in mold avoidance and yield stability along north-south rainfall gradients. In response to the high susceptibility of elite cultivars in subhumid zones of Senegal, we developed and characterized a recombinant inbred line (RIL) population derived from Nganda (grain mold-susceptible) and Grinkan (photoperiod-sensitive) varieties. The population was evaluated across three distinct agro-ecological zones over two years. Environmental indices derived from genotype-environmental interactions, together with defined growth windows, strongly influenced flag leaf appearance (FLA), a photoperiodic flowering trait. Plasticity parameters (intercept and slope) for environmental indices, FLA, grain mold severity, and yield enabled identification of loci contributing to flowering response, mold resistance, and yield stability. The maturity gene Ma1 and two QTLs for FLA, qFLA6.2 and qFLA6.3, were identified, stable across environments, and colocalized with grain mold and yield QTLs. The wild-type Ma1 allele from Grinkan delayed FLA and reduced grain mold damage but was not associated with increased yield. The Ma1 effect was confirmed using the developed breeder-friendly KASP marker, Sbv3.1_06_40312464K, in 174 F3 three-way cross families. Photoperiod-sensitive lines with intermediate-to-late FLA alleles showed strong negative associations with mold damage. Overall, the identified stable loci and candidate lines provide foundations for effective molecular breeding of climate-resilient varieties. PLAIN LANGUAGE SUMMARYGrain mold is a fungal disease that reduces sorghum grain yield and quality, particularly in subhumid climates. With the limited number of resistant elite varieties, photoperiod-sensitive flowering to day length variation can contribute to grain mold escape at the end of rainy seasons. We characterized 286 sorghum recombinant inbred lines across three contrasting environments over two years along rainfall gradients in Senegal. Using flag leaf appearance (FLA), which is a photoperiodic flowering trait, strong genotype-environment interactions for FLA and genotypic plasticity were revealed. We identified and validated the common genomic locus associated with FLA variation and its plasticity across environments, the canonical maturity gene Ma1, which was influenced by temperature variation across environments. The presence of Ma1 in the background of photoperiod-sensitive lines enhances grain mold avoidance and yield stability along rainfall gradients in Senegal. CORE IDEASO_LIWe investigated photoperiodic flowering plasticity in sorghum as a contributor to grain mold resistance and yield stability along rainfall gradients. C_LIO_LIThe Maturity locus Ma1 (qFLA6.1) is the major contributor of photoperiodic flowering and its plasticity across semi-arid and subhumid environments. C_LIO_LIHybrid genotypes carrying two stable loci qFLA6.1 and qFLA6.2 sustain high grain mold avoidance in diverse environments. C_LIO_LIPhotoperiod-sensitive lines with medium to late flowering times are effective in avoiding grain mold, while maintaining yield stability in subhumid regions. C_LI

Bioavailability of iron, an essential micronutrient to plants, is low in alkaline or calcareous soils, which are prevalent across semi-arid production regions. Breeding efforts to increase tolerance to iron deficiency chlorosis (IDC) in sorghum, a major crop of semi-arid regions, are confounded by spatial variation of stress severity in field trials. Here we developed and validated two high-throughput phenotyping approaches to address this challenge, with multi-spectral aerial imaging in the field and a controlled-environment assay to isolate the effects of iron bioavailability. In the field, severity and uniformity of stress are highly predictive of genetic signals for IDC tolerance (R2 > 0.6 for soil pH metrics and H2). Plot-level data filtering for stress conditions based on control genotypes successfully addresses field spatial variation (unfiltered H2 = 0.18 vs. filtered H2 = 0.4). The controlled-environment assay proxies field stress using iron sources with differential bioavailability, evidenced by high heritability ( H2 = 0.98) and phenotypic differential for hybrid control genotypes that matches field performance. Finally, we show that assay phenotypes are suitable for genome-wide association studies in global germplasm. Together, these field and lab phenomic approaches can be deployed to understand genetics of IDC tolerance and develop crops resilient to alkaline soils. HIGHLIGHTStress severity and uniformity greatly impact detection of genetic signals underlying iron deficiency chlorosis tolerance in sorghum. A controlled-environment assay reduces spatial heterogeneity and improves assessment of tolerance genetics.

DNA methylation plays a central role in the regulation of gene expression. In plants, methylation occurs in the CG, CHG and CHH contexts, via distinct DNA methyltransferases including MET1, CMT3 and the RNA-directed DNA Methylation (RdDM) pathway via DRM2. In interspecific hybrids, these epigenetic mechanisms are confronted to a mixed small RNA population and two subgenomes harbouring specific methylation patterns, therefore generating unique expression profiles. The aim of this work was to understand these regulations by analysing gene expression, DNA methylation and small RNAs in a Citrus hybrid resulting from the cross between C. reticulata (mandarin) and C. australasica (finger lime). Haplotype-resolved subgenomes assembly identified hundreds of allele-specifically expressed genes. Asymmetric reprogramming of methylation was observed, in particular an increase in CHH in C. australasica haplotype. Surprisingly, CHH methylation, usually associated with gene silencing, was correlated here with increased expression, but also 24nt small RNA populations at their promoter regions. Similar analyses of the parental lines and other citrus species suggest the correlation between CHH methylation-enriched promoter and high expression level is not due to the hybridization, but seem to be generally true for all citrus. These observations suggest that, in citrus fruit, RdDM could activate transcription. This work also provides a full pipeline to analyse the expression profiles and DNA methylation in complex hybrids, which could be crucial for anticipating varieties resistant to diseases and the current threats affecting citriculture such as the Huanglongbing disease.

Sesame (Sesamum indicum) is one of the earliest domesticated oilseed crops and is valued for antioxidant lignans that stabilize oil quality. However, the genomic and evolutionary history of the genus Sesamum, including the origin of its allotetraploid relative S. radiatum and the diversification of lignan metabolism, remains poorly understood owing to limited chromosome-scale genomic resources. Here we present chromosome-level genome assemblies for three wild Sesamum species, two Ceratotheca species and a Japanese sesame cultivar to reconstruct genome and karyotype evolution across the Sesamum-Ceratotheca complex. Comparative analyses show that the derived x=16 lineage originated from an ancestral x=13 karyotype through chromosome fission, fusion and translocation, whereas another x=13 lineage underwent extensive restructuring associated with retrotransposon expansion. Phylogenomics places Ceratotheca within the x=16 Sesamum clade and reveals that S. radiatum originated through hybridization involving a C. sesamoides-like ancestor. The antioxidative lignan gene CYP92B14 was reintroduced via the BB progenitor, linking hybridization with restoration of oil-stabilizing metabolism during sesame evolution.